Search for: All records

The tectonic setting of Northwestern South America is very complex due to interaction between the Nazca, Caribbean, and South America plates and several oceanic terranes. Of particular interest is the Nazca plate slab geometry change since the Mid-Miocene, and its tectonomagmatic effects on the overlying South American Plate. Several studies indicate that the modern Nazca slab is torn into two segments at 5.5oN in Colombia, the so-called "Caldas tear", where the northern segment dips at a shallow angle while the southern segment dips at a steep angle. This slab geometry difference is manifested in magmatic activity where only the region with steep slab shows typical subduction zone volcanism, while the one with flat slab lacks such activity. However, due to the limited number of geochronological data, the timing of this slab shallowing and tearing remains poorly understood. We conducted an extensive 40Ar/39Ar geochronology study on post mid-Miocene volcanic rocks in Colombia to study the temporal evolution of Nazca plate geometry and its effects on magmatism in Colombia. We find evidence of continuous magmatism north and south of the Caldas tear from 10.5 Ma to 6.4 Ma, with peak activity between 9-8 Ma. This is followed by a ~4 million year magmatic hiatus, until the resurgence of magmatism south of the Caldas tear as monogenetic domes at 2.1 Ma, and continuous volcanic activity in modern composite volcanoes since 1.1 Ma. Results of this study support the presence of a complex subduction system beneath Colombia where the northern segment slab has been flat since ~6.4 Ma ago, and the southern segment has re-steepened at ~2.1 Ma. This study showcases that the Nazca slab geometry was the most important factor driving magmatism in Colombia since the mid-Miocene. 

The tectonic setting of Northwestern South America is very complex due to interaction between the Nazca, Caribbean, and South America plates and several oceanic terranes. Of particular interest is the Nazca plate slab geometry change since the Mid-Miocene, and its tectonomagmatic effects on the overlying South American Plate. Several studies indicate that the modern Nazca slab is torn into two segments at 5.5oN in Colombia, the so-called "Caldas tear", where the northern segment dips at a shallow angle while the southern segment dips at a steep angle. This slab geometry difference is manifested in magmatic activity where only the region with steep slab shows typical subduction zone volcanism, while the one with flat slab lacks such activity. However, due to the limited number of geochronological data, the timing of this slab shallowing and tearing remains poorly understood. We conducted an extensive 40Ar/39Ar geochronology study on post mid-Miocene volcanic rocks in Colombia to study the temporal evolution of Nazca plate geometry and its effects on magmatism in Colombia. We find evidence of continuous magmatism north and south of the Caldas tear from 10.5 Ma to 6.4 Ma, with peak activity between 9-8 Ma. This is followed by a ~4 million year magmatic hiatus, until the resurgence of magmatism south of the Caldas tear as monogenetic domes at 2.1 Ma, and continuous volcanic activity in modern composite volcanoes since 1.1 Ma. Results of this study support the presence of a complex subduction system beneath Colombia where the northern segment slab has been flat since ~6.4 Ma ago, and the southern segment has re-steepened at ~2.1 Ma. This study showcases that the Nazca slab geometry was the most important factor driving magmatism in Colombia since the mid- Miocene. 

Abstract The iconic volcanoes of the Cascade arc stretch from Lassen Volcanic Center in northern California, through Oregon and Washington, to the Garibaldi Volcanic Belt in British Columbia. Recent studies have reviewed differences in the distribution and eruptive volumes of vents, as well as variations in geochemical compositions and heat flux along strike (amongst other characteristics). We investigate whether these along‐arc trends manifest as variations in magma storage conditions. We compile available constraints on magma storage depths from InSAR, geodetics, seismic inversions, and magnetotellurics for each major edifice and compare these to melt inclusion saturation pressures, pressures calculated using mineral‐only barometers, and constraints from experimental petrology. The availability of magma storage depth estimates varies greatly along the arc, with abundant geochemical and geophysical data available for some systems (e.g., Lassen Volcanic Center, Mount St. Helens) and very limited data available for other volcanoes, including many which are classified as “very high threat” by the USGS (e.g., Glacier Peak, Mount Baker, Mount Hood, Three Sisters). Acknowledging the limitations of data availability and the large uncertainties associated with certain methods, available data are indicative of magma storage within the upper 15 km of the crust (∼2 ± 2 kbar) beneath the main edifices. These findings are consistent with previous work recognizing barometric estimates cluster within the upper crust in many arcs worldwide. There are no clear offsets in magma storage between arc segments that are in extension, transtension or compression, although substantially more petrological work is needed for fine scale evaluation of storage pressures. 

The chemistry of erupted clinopyroxene crystals (±equilibrium liquids) have been widely used to deduce the pressures and temperatures of magma storage in volcanic arcs. However, the large number of different equations parameterizing the relationship between mineral and melt compositions and intensive variables such as pressure and temperature yield vastly different results, with implications for our interpretation of magma storage conditions. We use a new test dataset composed of the average Clinopyroxene-Liquid (Cpx-Liq) compositions from N = 543 variably hydrous experiments at crustal conditions (1 bar to 17 kbar) to assess the performance of different thermobarometers and identify the most accurate and precise expressions for application to subduction zone magmas. First, we assess different equilibrium tests, finding that comparing the measured and predicted Enstatite-Ferrosillite and KD (using Fet in both phases) are the most useful tests in arc magmas, whereas CaTs, CaTi and Jd tests have limited utility. We then apply further quality filters based on cation sums (3.95–4.05), number of analyses (N > 5) and the presence of reported H2O data in the quenched experimental glass (hereafter ‘liquid’) to obtain a filtered dataset (N = 214). We use this filtered dataset to compare calculated versus experimental pressures and temperatures for different combinations of thermobarometers. A number of Cpx-Liq thermometers perform very well when liquid H2O contents are known, although the Cpx composition contributes little to the calculated temperature relative to the liquid composition. Most Cpx-only thermometers perform very badly, greatly overestimating temperatures for hydrous experiments. These two findings demonstrate that the Cpx chemistry alone holds very little temperature information in hydrous systems. Most Cpx-Liq and Cpx-only barometers show similar performance to one another (mostly yielding root mean square errors [RMSEs] of 2–3.5 kbar), although the best Cpx-only barometers currently outperform the best Cpx-Liq barometers. We also assess the sensitivity of different equations to melt H2O contents, which are poorly constrained in many natural systems. Overall, this work demonstrates that Cpx-based barometry on individual Cpx only provides sufficient resolution to distinguish broad storage regions in continental arcs (e.g. upper, mid, lower crust). Significant averaging of Cpx compositions from experiments reported at similar pressures can reduce RMSEs to ~1.3–1.9 kbar. We hope our findings motivate the substantial amount of experimental and analytical work that is required to obtain precise and accurate estimates of magma storage depths from Cpx ± Liq equilibrium in volcanic arcs. 

Understanding the processes that initiate volcanic eruptions after periods of quiescence are of paramount importance to interpreting volcano monitoring signals and mitigating volcanic hazards. However, studies of eruption initiation mechanisms are rarely systematically applied to high-risk volcanoes. Studies of erupted materials provide important insight into eruption initiation, as they provide direct insight into the physical and chemical changes that occur in magma reservoirs prior to eruptions, but are also often underutilized. Petrologic and geochemical studies can also constrain the timing of processes involved in eruption initiation, and the time that might be expected to elapse between remote detection of increased activity and eventual eruption. A compilation and analysis of literature data suggests that there are statistical differences in the composition, volume, style and timescales between eruptions initiated by different mechanisms. Knowledge of the processes that initiate eruptions at a given volcano may thus have significant predictive power. 

Luo et al.1 present exciting new data on lunar basalt samples erupted at ~2 Ga, and brought to Earth by the Chang’e-5 (CE-5) mission. These samples offer important new opportunities to understand lunar magmatic systems. Luo et al.1 use Clinopyroxene-Liquid (Cpx-Liq) thermobarometry and pMELTS modelling of mineral compositions to determine the pressures (P) and temperatures (T) of magma storage on the moon. However, in this comment we discus the large analytical errors associated with their measurements, and the implications for their intepretation. 

In 1918, the geologist Emile Grosse was commissioned to conduct geological studies in the Amagá Basin, Antioquia, Colombia. In 1923, Grosse finished a comprehensive cartographic work that became the cornerstone for the geology of the northwest (NW) Colombian Andes. Today, 100 years later, the volcanoclastic strata preserved in the Amagá Basin are crucial for understanding major Oligocene to Pliocene tectonic events that occurred in the NW South-American margin, including the fragmentation of the Nazca Plate, the collision of the Panamá-Chocó Block, and the shallowing of the subducted slab. Our contribution includes new mineral chemistry and zircon petrochronological data from the Combia Volcanic Complex and published data to provide a review of the Oligocene to Pliocene deformation, sedimentation, and magmatic patterns in the Amagá Basin and their implications for the tectonic evolution of NW South America. The Amagá Basin was the result of the Eocene to Oligocene uplift of the Western Cordillera followed by the Middle Miocene to Pliocene uplift of both the Central and Western cordilleras, events that modified the Miocene drainage network in the Northern Andes. Coeval with the final Miocene deformation phases in the Amagá basin, the magmatism of the Combia Complex was the result of subduction magmas emplaced in a continental crust affected by strike-slip tectonics. 

Abstract The composition of clinopyroxene and clinopyroxene-liquid (Cpx-Liq) pairs are frequently used to calculate crystallization/equilibration pressures in igneous systems. While canonical uncertainties are often assigned to calculated pressures based on fits to calibration or test datasets, the sources of these uncertainties (and thus ways to reduce them) have not been rigorously assessed. We show that considerable uncertainties in calculated pressures arise from analytical error associated with Electron Probe Microanalyser (EPMA) measurements of Cpx. Specifically, low X-ray counts during analysis of elements with concentrations <1 wt% resulting from insufficient count times and/or low beam currents yield highly imprecise measurements (1σ errors of 10–40% for Na2O). Low analytical precision propagates into the calculation of pressure-sensitive mineral components such as jadeite. Using Monte Carlo approaches, we demonstrate that elemental variation resulting from analytical precision alone generates pressures spanning ~4 kbar (~15 km) for a single Cpx and ~6 kbar for a single Cpx-Liq pair using popular barometry expressions. In addition, analytical uncertainties in mineral compositions produce highly correlated arrays between pressure and temperature that have been previously attributed to transcrustal magma storage. Before invoking such geological interpretations, a more mundane origin from analytical imprecision must be ruled out. Most importantly, low analytical precision does not just affect the application of barometers to natural systems; it has also affected characterization of Cpx in experimental products used to calibrate and test barometers. The impact of poor precision on each individual measurement is often magnified by the small number of measurements made within experimental charges, meaning that low analytical precision and true variability in mineral compositions have not been sufficiently mediated by averaging multiple EPMA analyses. We compile the number of Cpx measurements performed in N = 307 experiments used to calibrate existing barometers, and N = 490 new experiments, finding ~45% of experiment charges were characterized by ≤5 individual Cpx analyses. Insufficient characterization of the true composition of experimental phases likely accounts for the fact that all Cpx-based barometers exhibit large errors (± 3 kbar) when tested using global experimental datasets. We suggest specific changes to analytical and experimental protocols, such as increased count times and/or higher beam currents when measuring low concentration elements in relatively beam resistant Cpx in experiments and natural samples. We also advocate for increasing the number of analyses per experimental charge, resolving interlaboratory analytical offsets and improving data reporting. Implementing these changes is essential to produce a more robust dataset to calibrate and test the next generation of more precise and accurate Cpx-based barometers. In turn, this will enable more rigorous investigation of magma storage geometries in a variety of tectonic settings (e.g. distinguishing true transcrustal storage vs. storage in discrete reservoirs). 

We present Thermobar, a new open-source Python3 package for calculating pressures, temperatures, and melt compositions from mineral and mineral-melt equilibrium. Thermobar allows users to perform calculations with >100 popular parametrizations involving liquid, olivine-liquid, olivine-spinel, pyroxene only, pyroxene-liquid, two pyroxene, feldspar-liquid, two feldspar, amphibole only, amphibole-liquid, and garnet equilibria. Thermobar is the first open-source tool which can match up all possible pairs of phases from a given region, and apply various equilibrium tests to identify pairs from which to calculate pressures and temperatures (e.g. pyroxene-liquid, two pyroxene, feldspar-liquid, two feldspar, amphibole-liquid). Thermobar also contains functions allowing users to propagate analytical errors using Monte-Carlo methods, convert pressures to depths using different crustal density profiles, plot mineral classification and mineral-melt equilibrium diagrams, calculate liquid viscosities, and convert between oxygen fugacity values, buffer positions and Fe speciation in a silicate melt. Thermobar can be downloaded using pip and extensive documentation is available at https://thermobar.readthedocs.io/.